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Abstract 

It was recently shown that in cosmology the gravitational action of 
faraway matter has quite relevant effects, if retardation of the forces 
and discreteness of matter (with its spatial correlation) are taken into 
account. Indeed, far matter was found to exert, on a test particle, 
a force per unit mass of the order of 0.2 cHq. It is shown here that 
such a force can account for the observed rotational velocity curves 
in spiral galaxies, if the force is assumed to be decorrelated beyond a 
sufficiently large distance, of the order of 1 kpc. In particular we fit 
the rotation curves of the galaxies NGC 3198, NGC 2403, UGC 2885 
and NGC 4725 without any need of introducing dark matter at all. 
Two cases of galaxies presenting faster than keplerian decay are also 
considered. 

Keywords: galaxies: general - galaxies: kinematics and dynamics - 
galaxies: fundamental parameter 

1 Introduction 

Perhaps, the actual problem of gravitation theory is to justify the large speeds 
observed in objects of galactic or larger sizes. This problem became increas- 
ingly puzzling, starting from the end of the years seventies. Indeed precise 
measurements (with a variety of techniques) of the speed of gases in spiral 
galaxies made apparent that the action of the galaxy on the gas, computed 
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according to the standard theory of gravitation, cannot account for the mea- 
sured velocities, if the gas is assumed to be a stable component of the galaxy. 
In the literature, while a lot of efforts are devoted to computing the grav- 
itational action of a galaxy on the nearby gas, no attention is paid to the 
gravitational effects exerted by the other galaxies, particularly the far away 
ones. This attitude, of neglecting the faraway matter, would indeed be cor- 
rect if the other galaxies were uniformly distributed. However, up to the 
present limit of observations (at least 300 Mpc, according to Gabrielli et al. 
|1999[ 12004] : Sylos Labini et al. [1998]), it was suggested that the galaxies 
may instead be distributed according to a very complex fractal structure. It 
was shown by Carati et al. ( |2008j ) that the gravitational action of faraway 
matter computed according to the General Relativity has quite relevant ef- 
fects, if both the retardation of the forces and the discreteness of matter 
(with its spatial correlation) are taken into account. Indeed, the far matter 
was found to exert, on a test particle, a force per unit mass of the order 
of 0.2 cHq, where Hq is the Hubble constant, and c the speed of light. It 
was also shown how the gravitational force due to the faraway matter may 
account for the speed of the galaxies in clusters. 

The aim of the present paper is to estimate the effect that the faraway 
matter produces on the rotation curves in spiral galaxies. We show that such 
an effect becomes very relevant if one assumes that the gravitational field 
due to the faraway matter is decorrelated beyond a sufficiently large distance 
(of the order of 1 kpc). So such a force goes unnoticed for objects of the solar 
system's size (or smaller), while starts becoming important from objects of 
galactic size on, being the predominant one for clusters of galaxies and larger 
structures. In particular, we show that the rotation curves of four galaxies 
(with very different masses and sizes) can be fitted taking into account such 
a force, leading to reasonable values of the luminosity-mass ratio (see Table 
I of Section 4). Two galaxies presenting a faster than keplerian decay are 
also considered. 

The paper is organized as follows. In Section 2, we recall the model 
that was introduced by Carati et al. ( |2008] ) in order to deal with the 
gravitational effect of faraway matter. In Section 3 the effect on the rotation 
curves is evaluated, and in Section 4 the fit to rotational curves for the 
galaxies NGC 3198, NGC 2403, NGC 2885, NGC 4725 is reported. Two 
cases of rotation curves presenting a faster than keplerian dacay are also 
briefly discussed. Two Appendices, in which some technical computations 
are reported, complete the paper. 



2 



2 Definition of the model 



We briefly recall in this section the model that was introduced by Carati et 
al. ( j2008] ) in order to estimate the gravitational action of far away matter. 
To this end, flrst of all one has to estimate the gravitational field generated 
by the distant galaxies (thought of as point masses) considered as sources 
in Einstein's equation, in the approximation in which their positions and 
velocities are assigned, in the simplest way compatible with observational 
cosmology. 

To this end, the positions of the galaxies are considered as random 
variables, whose statistical properties will be discussed later. The velocity 
field of the sources is instead taken according to Bubble's law, the peculiar 
velocities being altogether neglected. Namely, taking a locally Minkowskian 
coordinate system centered about an arbitrary point, a particle with position 
vector q is assumed to have a velocity 



where, for the sake of simplicity of the model, the Hubble constant Hq is 
assumed to be time independent. It is easily established that the considered 
chart has a local Hubble horizon _Ro = c/Hq, where the galaxies' speed equals 
that of light. 

The energy-momentum tensor T'^'^ is then given by 



where is the number of galaxies, Mj and 7^ are the mass and the Lorentz 
factor of the j-th one, g is the determinant of the metric tensor (which is 
considered as an unknown of the problem), 5 the Dirac delta function, and 
the dot denotes derivative with respect to proper time along the worldline of 
the source. 

The study of the solutions of Einstein's equations with the energy-momentum 
tensor ([2]) as a source is a formidable task, and so the study of Carati et al. 
( |20U8] ) was performed within the limits of a perturbation approach, consid- 
ering the energy-momentum tensor T^'^ (|2]) as a perturbation of the vacuum. 
Following the standard procedure (see Einstein |1922] or Weinberg |1972] ). 
one then has to determine a zero-th order solution (the vacuum solution), 
and solve the Einstein equations, linearized about it. The simplest consistent 
zero-th order solution is the fiat metric because, coherently, the perturbation 
is then shown to be small (at least if the parameters of the model, such as 
the density, are chosen in agreement with the observations). 



q = ^^0 q , 



(1) 
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Thus the metric tensor g^^, is written as a perturbation of the Minkowskian 
background r/^^, namely, as g^i^ = rj^^ + h^^, and it is well known that in the 
linear approximation the perturbation h^y then has to satisfy essentially the 
wave equation with T^^^ as a source. More precisely, one gets 



where G is the gravitational constant, h the trace of /i^j^, and □ = {1/0^)8^ — 
A2. The solutions are the well known retarded potentials 

(4) 



c4 ^ 7,- |x-q,-| 



t=Uct 



(with q'^-'^ = qj). 

In the present case in which we are concerned with the rotation curves 
of spiral galaxies, one has also to take into account the contribution of local 
matter, so that, in the spirit of the considered approximation, the total metric 
tensor is given by the sum of that due to the local matter, and of that due 
to the distant one. For the local metric one furthermore makes the 
Newtonian approximation, i.e., all components are assumed to vanish, apart 
from h^QQi which is set equal to the Newtonian potential V^°'^{yi) /2c^ due to 
local masses. By local matter we mean the one contained in the galaxy, so 
that we forget, in the rest of the paper, the possible contribution of the dark 
matter halo. The total metric is thus written as 

2y'°^(x) 

g^lu — Vfj^i^ "I Ooi/Oo^t + iT'fii^ ) 

with h^i, given by (jl]). 

So the equation of motion for a test particle contains, in addition to 
the Newton force produced by the mass inside the galaxy, coming from the 
term h^°^, also some terms coming from the field due to the distant matter. 
How large is this additional term with respect to the local Newton force? 
The answer depends, on the assumptions made on the positions q^- of the 
constituents of the far away matter which appear in formula (jl]). 

Carati et al. ([2008]) assumed the positions q^- to be random variables, 
having a probability distribution which is isotropic (rotationally invariant). 
A simple computation then shows that the mean force (per unit mass) acting 
on a test particle vanishes. Thus the order of magnitude of the force is given 
by its standard deviation. 

The estimate of the standard deviation requires an additional hypotesis on 
the distribution of the quantities q^-. If they are assumed to be independent 
and identically distributed, then the force turns out to be negligible. 
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On the other hand, it is known that the positions of different galaxies are 
correlated, having a distribution with a fractal character, and it has been 
suggested that this may happen up a certain quite large distance (see Sylos 
Labini et al. [T998] : GabrieUi et al. plKMl [1199]). Carati et al. f plKlH] ). 
assuming that the distribution is fractal at all distance scales, show that the 
force per unit mass due to the distant galaxies is of the order of 0.2 cHq. 

Now, as pointed out by Milgrom ( |1983[ 11987] ). it is precisely when the 
acceleration of the stars is of the order of cHq that the rotation curves of 
the spiral galaxies depart from the behaviour expected from the Newtonian 
force due to local visible matter. Our purpose is then to compute the effect 
of the far field, precisely in such a region. It will be shown that such an 
the effect can be described as corresponding to an effective potential which 
acts on the stars by modifying the profile of the rotation curves as to fit the 
actual observed one. 

3 The effect of far away matter on the rota- 
tion curves 

In this Section, we estimate the effect of the field produced by the far away 
matter on the rotation curves. So, first of all, one has to estimate the influence 
of the far field on the motion of stars. 

In general, as there is no reason to expect that the far field should be 
spherically symmetric, its first effect amounts to destroy the circular orbits, 
from which every discussion about the rotation curves usually starts. To 
tackle this poblem, one can imagine that the effect of the far field is not too 
large, and that the orbits are still close to circular ones. In this case one can 
average the radial component of the equation of motion along a circle, thus 
obtaining an effective radial equation from which the angular velocity of the 
stars as a function of the radius may be estimated. 

Using cylindrical coordinates (r, 6, z), and considering only motions which 
lie in the galactic plane z = {so that z =z = 0), the radial equation of 
motion turns out to be 



In particular, in terms of the cartesian components /i^i/ given by (jl]) and of 




(5) 
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the local potential h^°^, the relevant components are given by the formulas 



9rr 

9rd 
90r 



9oe 
900 



(l + hii cos^ 9 + /i22 sin^ 9 + hi2 sin 9 cos 6') 

'^((^ii ~ ^22) sin 6* cos 6* + hi2 cos 26*) =^ rV 
hoi cos 6* + ho2 sin 6^ 

2 (1 + hu sin^ ^ + h22 cos^ ^ + sin ^ cos = r'^{l + A) 



r 



r 



(/io2 cos 9 — hoi sin 9) =^ ri3 



(6) 



The quantities A, B and "D, implicitly defined in such formulas, will be used 
later. In our setting, equation ([5]) is a stochastic one, because of the stochastic 
character of the terms /i^j,, and it is very hard to be dealt with. 

To reduce our problem to a tractable one, we make essentially two hy- 
potheses. First of all we suppose that it is meaningful to average the equation 
with respect to 6*, i.e., we suppose that 6^ is a fast variable with respect to r, 
and consequently the derivative r must be thought of as "small". Thus, in 
the averaged equation, we will forget all terms proportional to r. 

Now, as the galaxies do not change their dimensions on the time scale 
of the revolution of a typical star, one can readily give an estimate of the 
term d{grrr)/dt of ([S]), which turns out to essentially vanish. Obviously, also 
the terms which are the derivatives with respect to 9 of components of the 
metric tensor will have a vanishing averages. One is therefore left with the 
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following equation in the unknown v = r9: 



+ 



27rr 



{A + rdrA) dl 



V 

r 



— Lb + rdrB) dl- — (I) dtV dl 
.2'Kr I 27rr 



V 

J r 



2'nr 



dt9or dl + — (f) drflQQ dl 



(7) 



where use was made of expressions ([6]) for the components of the metric 
tensor. This equation has to be compared (see for example Bertin [2UUUj ) 
with the familiar one 

which is obtained when the only acting force is assumed to be the Newtonian 
one due to local matter. So, at first sight, the structure of the equation seems 
to have changed just a little bit: there appears a term which is linear in f , 
but this modification turns out to be negligible. However, a more relevant 
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contribution (in some hypotheses to be discussed in a moment) is given by 
the term 

^ jdtg,r^l = -drV'ff . (8) 

In this respect, the faraway field acts as to produce a modification of the 
local potential, giving rise to a sort of effective potential; this effect might 
also be mimicked by a local distribution of matter (of an unknown type). 

The relevant difference between our formula ([7]) and the standard one, 
consists in the fact that in ([7]) the terms depending on the far fields have a 
stochastic character, so that, as already pointed out, we cannot predict their 
exact values, but can only estimate their order of magnitude by looking at 
their average and variance. 

Notice that, as we are integrating along a circular path, the magnitude 
of the random terms will depend on the correlation of the field at different 
points of the path. Being interested only in the variance, it will suffice to 
just consider the two-point correlations. 

In principle such correlations should be computed from the expression (jlj) , 
but at the moment we are unable to do this in a consistent way. Therefore we 
make here our second assumption, i.e., it will be assumed that the correlations 
decay in an exponential way on a certain lenght scale /, namely, we assume 

{dagt,v{x)d„g^y{y)) = exp - ^^-y^) ( (9^^^z.(a;))^) . (9) 

Here / is considered as a free parameter, to be fitted by comparison with 
the actual observations of galaxy rotation curves. If one obtains consistent 
values of / for different galaxies, the ansatz ([9]) might be considered as sound. 
If / is large, the field can be considered as smooth, becoming more stochastic 
as / decreases. We will find for I a value of 1 kpc, so that the far field can be 
considered as stochastic for what concerns the motion of stars in a galaxy, 
and smooth instead, for example for what concerns the motion of planets 
around a star. 

We recall that, via the Wiener-Khinchin formula (see Appendix B), one 
can obtain the correlation of a function, from the correlations of its deriva- 
tives. In our case, we obtain 

(^7..(x)M2/)> = ^(l-e-^(l + ^))((5.M^))'> • (10) 

Using ([9]) and f lTO|) . the averages and variances of all terms involved in equa- 
tion ([7]) can be computed as shown in Appendix A. The results of such 
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computations are the following. Denote 



A 



C 



def 



def 



def 



-^IaAI, A2 = T^(fdrAdl, 
Inr J ZTrr J 

Zirr J Zirr J 

— I dtgor dl , = — / dtV dl . 

2'Kr J ' 27rr/ 



Then, for the averages and the variances one finds the following results: 



where the constants are given by 7a — 0.05, 7^ ~ 0.026 and 7D ~ 0.08. 

Notice that in the first line the variance is not computed because, being 
the average not vanishing, it is just the average which gives the order of 
magnitude of the terms in equation ([7]). Using relations f lTTj) . for r in the 
range of galactic distances, i.e., r < 30 kpc, one checks that all terms in 
([7]) are negligible, except for three of them: the local Newtonian potential, 
the quadratic term in the tangential velocity v, and the term containing 
f dtgordl/2'Kr. Thus takes the simpler form 



Here, the value indicated after the sign ~ is just the standard deviation. In- 
deed, as explained in the Introduction, we are dealing with random variables 
having zero mean, so that their order of magnitude is given by their standard 
deviation. Remember that the sign can be positive or negative: if the term 
is positive it acts as a pressure and so it helps keeping the galaxy stable, 
while the opposite occurs if it is negative. So one can conjecture that the 
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Figure 1: Plot of the rotational velocity v vrtsus distance r from the galactic 
center, for the galaxy NGC 3198. Squares (□) are the velocities determined 
from HI observations (see Bergman }1989j ). solid line is the theoretical curve 
(see text). Dashed line is the contribution due to the local matter. 



9 



positive case occurs more frequently in observations, but it would be possible 
in principle to observe also cases in which the term has a negative value. 

We show in the next section that, if the term is positive, the magnitude 
of the correction due to the far field is able to flatten the rotation curve in 
the external region of the galaxy. 

4 Fitting some observed rotation curves 

To make a comparison between the predictions given by (|T2|) . and the ob- 
served rotation curves of galaxies, one should have available an expression for 
V''°'^{r), or, equivalently, one should assign the local distribution of matter. 
In the literature, there are several papers which illustrate different ways in 
which the distribution of matter can be estimated, starting for example from 
the measured distribution of luminosity of the galaxies (see Toomre |1963] : 
Freeman [1970] : Carignan&:Freeman |1985j : Carignan et al. |1988j ). 

Table 1: Values of total mass, mass-luminosity ratio and correlation length 
for the four galaxies considered. 



Galaxy 


Mass 


Mass/Luminosity 


Correlation length 


NGC 3198 


4.0 10^° Mq 


4.6 Mq/Lq 


0.6 Kpc 


NGC 2403 


3.5 10^° Mq 


4.4 Mq/Lq 


0.8 Kpc 


NGC 2885 


1.0 10^2 Mq 


2.1 Mq/Lq 


1.7 Kpc 


NGC 4725 


1.1 10^1 Mq 


2.1 Mq/Lq 


3.1 Kpc 



We take instead the simpler path which consists in assuming a functional 
form for V''°'^{r), with a minimal number of parameters (essentialy two, total 
mass and radius of the galaxy) and then trying to determine these parameters 
by a best fit with the observed rotation curves. In particular we concentrate 
on the rotation curves of the four galaxies NGC 3198, NGC 2403, UGC 2885 
and NGC 4725. 

For what concerns the local potential V'°'^(r), in the literature one often 
makes use of potentials of the form 



first introduced by Miyamoto&Nagai ([ 1975] ). Here M is the total mass of 
the object, to be understood in the following sense. First one introduces the 
matter density p from the Poisson equation A2V''°^ = AnGp, then one easily 
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Figure 2: Plot of the rotational velocity v versus distance r from the galactic 
center, for the galaxy NGC 2403. Squares (□) are the velocities determined 
from HI observations (see Bergman |1987] ). while triangles (A) are Ho; ob- 
servations (see Blais-Ouellette et al. |2004] ) . The solid line is the theoretical 
curve (see text), while the dashed line is the contribution due to the local 
matter. 
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Figure 3: Plot of the rotational velocity v versus distance r from the galactic 
center, for the galaxy NGC 4725. Squares (□) are the velocities determined 
from HI observations (see Wevers et al. |1984j ). The solid line is the theo- 
retical curve (see text), while the dashed line is the contribution due to the 
local matter. 
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Figure 4: Plot of the rotational velocity v versus distance r from the galactic 
center, for the galaxy UGC 2885. Squares (□) are the velocities determined 
from HI observations (see Roelfsema& Allen |1985] ). while triangles (A) are 
Ha observations (see Rubin et al. [1985] ). The solid line is the theoretical 
curve (see text), while the dashed line is the contribution due to the local 
matter. 
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checks that the integral of p, over the whole space, is equal to the parameter 



For what concerns the other two parameters a and b, one has first of 
all that the ratio b/a determines the flatness of the mass distribution p: for 
b/a > 1, the distribution is essentially spherical, while it reduces to a singular 
disk if the ratio vanishes. The parameter a (at fixed ratio b/a) is then the 
lenght scale of the distribution p. 

So, for the galaxies NGC 3198 and NGC 4725 we take the potential 
yioc^^^ in the the Miyamoto-Nagai form, supposing that the only relevant 
contribution to the potential comes from the star disk. In the case of the 
other two galaxies, i.e, NGC 2403 and UGC 2885, it appears that there is also 
a non negligible contribution from the inner part of the galaxy (the so called 
bulge), so that we take the local potential as the sum of two Miyamoto-Nagay 
potentials, with different parameters, i.e., we take 



the first term referring to the bulge, the second to the disk. In this case, the 
total mass M of the galaxy will be given by M = Mi + M2. 

Now, using the above expressions for the local potentials, we made a best 
fit to the actual observed curve in order to determine both the parameters 
of the potential and the correlation length I. The resulting rotation curves 
are reported as solid lines in Figures [IHH The dashed lines are instead the 
rotation curves which would be obtained if the contribution of the distant 
matter were neglected. As one sees, the contribution of far matter is essential 
to keep fiat the rotation curves at large distances. The values of the masses 
M (in solar units Mq), the mass luminosity ratios (in solar mass over solar 
luminosity Lq units), and the correlation lengths / in kiloparsecs for the four 
galaxies are reported in Table [1] 

At variance with the effect due to dark matter, the far away matter could 
in principle also produce a negative pressure, because, as already explained, 
the term d^V'^^^ in equation f lT^ can be positive or negative. This would 
imply that the rotation curve decays, beyond the galaxy's edge, faster than 
keplerian. By a quick search, we found in literature at least two cases of 
this behavior: the rotation curve of the galaxy NGC 864, (see Espada et al. 
[2005] ) and that of galaxy AGC 400848 (see Catinella et al. [2005] ). We are 
able to fit also such curves, as shown in Figures O and El with again a value 
of the correlation lenght / of the order of 1 kpc. 

In any case, these curves are shown here only for the sake of illustration: 
in fact, first of all the contribution of the local matter was computed with 
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Figure 5: Plot of the rotational velocity v versus distance r from the galactic 
center, for the galaxy NGC 864. Squares (□) are the velocities determined 
from HI observations (see Espada et al. |2005j ). solid line is the theoretical 
curve (see text). 
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Figure 6: Plot of the rotational velocity v versus distance r from the galactic 
center, for the galaxy AGC 400848. Squares (□) are the velocities determined 
from Ha observations (see Catinella et al. [2005] ). solid line is the theoretical 
curve (see text). 
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an "ad hoc" potential without a deep confrontation with the data taken 
from luminosity to see if they are compatible. More importantly, we did not 
investigate whether there are some perturbations coming from known objects 
(such as nearby galaxies), which could perturb the velocity profile. 

5 Conclusions 

It was shown here that the faraway matter can exert a force of such a magni- 
tude as to account for the observed rotation velocity curves in spiral galaxies, 
if the force is assumed to be decorrelated beyond a sufficiently large distance, 
of the order of 1 kpc. In particular we fit the rotation curves of the galax- 
ies NGC 3198, NGC 2403, UGC 2885 and NGC 4725 without any need to 
introduce dark matter at all. 

We also showed that in principle faraway matter could act as a negative 
pressure, which would steepen the fall of of the rotation curves (with respect 
to the keplerian decay), as apparently observed for some galaxies. We con- 
sidered, in a preliminary way, two such cases, for which our approach seems 
to apply, giving a quite good fit. This subject is however left for a future 
work. 
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A Computation of variances 

In this appendix we compute the average and the variance of the terms which 
appear in the equation of motion ([7]). We recall that, from the fourth of (jH]) 
of Section 3 and from the expression of h^y given bylHof Section 2, one has 

A_ 4^^,^ (#cosg + #sing)' 2G ^ M, 
^~ ^^^^^^ \^Jy^\ V^|qO-)-x| ' ^^^^ 

where we are using cylindrical coordinates (r, 6, z), and the point x at which 
we are evaluating the field lies on the galactic plane z = 0. 

We have now to consider -Ai = ^ j Adl, and in particular (^i). How- 
ever, one obviously has 

Inr 
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because (A) is independent of the angle 9, due to the isotropy of the q^-^^ 
distribution. So we can take 6 = 0, and we obtain 

i j 



^4 v^-^J^-^J— ^J-- -Ji ^2 ^^|q(i)| 



where Rj = is the distance of the source q'--'^ from the galaxy center, 

9j is the angle between the source q*^-^^ and the galactic plane, and (pj a third 
angle on the galactic plane. We have obviously omitted any term of order 
|x|/i?j. To compute the average, we can take the continuum limit (as was 
done in paper Carati et al. [2008]), setting Mj = peffdV, so that one gets 

( A ) = -^^^(2H^, J R'sm''<Psm^9dRded^ + 

c'J RsmedRded<p) = -^p,^^ffX(f + ^) = ~l , 



where we have used the expression (8) of Carati et al. ( |2008] ) for the effective 
density Pe//- This gives the first of ( ITT]) . 

We come now to the second of (11 II) . In the first place one has to consider 
that, in computing the derivative which appears in the forces, one gets many 
terms, the largest of which comes from the fact the time at which we take 
the position of the sources is retarded, i.e. is a function of x. One has then 

46" , 2 (qi^^ cos e + 4^'^ sin 6) (qi^^ cos 6 + q^^ sin 9) 

drA = ^> ' I I 

^ ^ qO) - X 

i 

2Gv- ^7 ^Jq^-^'^-xl 

> — -7TT h smaller terms , 



c — ' q 

3 ' 



where "smaller terms" has to be understood in the sense that they give a 
smaller contribution to the mean and to the variance. Now one has 

(gi-'^ cos6' + gj^^ sin 61) 



a.|q(^)-x| = g.(g(^)+r ^^- ^"^^ ""'^ 0((r/g(^-^)^)) 

+ 0(r/g(^')) , 



{0^ cos Q + q;^^ sin 



|qO-)| 
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where we have used the fact that q^^^ ^ r. We get then, using Hubble's law, 
the following expression 

9^^_ggg!^M/(^"'-;,tf-'°'')°. (14) 
c . q 

where, from now on, we consider only the most relevant terms in the ex- 
pressions. Now, due to the rotational invariance of the distribution, one gets 
immediately 

(-^ (fdrAdl) = l{drA)dl = , 



*27rr / 27rr _ 

i.e. the first equality of the second line in f ITT]) . The estimate of the variance 
(Tyi2 is obtained in the following way. By definition one has 

^-42 = (^^^y d^i^" dl2drA{x)drA{y) 
and using the hypothesis on the correlations one gets 



Now, due to the rotational invariance of the distribution, the average inside 
the integral does not depend on the angle 6 and can be taken out from the 
integral. So one gets 



= f dh I dOexp ( - - sine/2) = 7i-((9.A) 



Now, the factor 71, which is essentially the value of the inner integral, depends 
indeed on the ratio r/l, but in the range of interest is essentially equal to its 
asymptotic value l/vr (such a value can be computed for example by applying 
Laplace method to the inner integral). So we simply put 



Z"^'" / 2r \ 1 / 

/ d^exp -—sine/2) = — 
Jo V i / 71 r 



having stipulated that Z/r is sufficiently large. 

We are left with the computation of the variance {{drA)"^), which can be 
performed as follows. From relation f|T^ . taking the square and using the 
rotational invariance of the distribution to compute it for ^ = 0, one gets the 
following expression 



j,k 



2 
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with Mtot = NM, where N the total number of galaxies, while 7^ is a 
numerical coefficient which depends on the distribution of the galaxies. If 
we choose the distribution as discussed in Section 2, one can estimate it by 
numerical computations, which give 7^ ~ 0.05. Now, using for Mtot the 
expression 

Mtot = iA/3)7Tp,ffRl , (15) 
with peff as given by (8) of Carati et al. ( [2008j ). one eventually finds 

^aI fHoV 

(^A2 = 

vr r \ c / 

which is the second relation in the second line of fllip . 

Let us now come to the third and fourth lines of (|TT|) . The average of Bi 
vanishes by symmetry, as can be seen from its expression 

c3 27rr/^ ^ |qO)| 

where use was made of Bubble's law. For what concerns the variance 
one has 

= ^ / / d^2(i5(x)i3(y)) 



27r 



/ ^ (^1 - e-^T f (1 + ; sin -) ) 

Jo 2rsine/2 V ^ I 2V 



where we used f[T9|) of Appendix B to estimate the correlation of a function 
in terms of the correlations of its derivatives. Use was also made of the 
isotropy of the distribution, so that (VS^)) turns out to be constant along 
the integration path. The integral in the second line can be computed nu- 
merically, and with a good accuracy (in the range of interest) one finds the 
value 27r/^(Z/r)°-^ so that one has 

aB^^{{yBf)l^(^'' . (16) 

Exactly in the same way one can check that the average of B2 vanishes, while 
the variance is given by 

^ {[d^Bf)- . (17) 
Trr 
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The expression of d^B, retaining only the largest terms, is 

,(i)|2 



_4G ^ (gj^'^ cos ^ - gi^'^ sin 9) jgi'^ cos g + g^^ sin g)) 



so that, using Hubble's law, one gets 



^^^^^^^-\-^) ^1^ |qO-)|2|qW|2 7^ ) 



Here, as before, one has Mtot = NM where N is the number of galaxies, while 
the numerical coefficient 7s depends on the distribution of the galaxies; for 
our fractal model a numerical estimate gives 7^ — 0.026. Again, inserting 
this expression into relation ( IT7|) . and using for Mtot the value ( 1TB]) . one gets 
the fourth line of relation f llip . The third line is obtained simply by observing 
that 

((VSf)=3((9,.i3)^) . 

We come now to the fifth line of (ITT]) . One has 

4(7 v-^ , ^ q^J^ cos e + gj^ sin 6 
^9r0 = 2^ —y-— , 

j 

so that one gets (forgetting as usual all small terms) the relation 
^ AG Hi ^ ^ ^ gj^'^ cos 9 - sin Q 

j 

c2 ^ ^ IqO') - x| 

i 

From this expression it follows immediately that (dtgro) = so that one has 
also 



(C) = -^/(c%,o)d/ = 
Znr J 



which is the first relation of the last line in fllip . The computation of the 
variance cx^ is made in exactly the same way as was done before for all other 
variances, obtaining 

^2 



= j j dl2{dtgrf){^)dtgrM) 



2tx 



/ d^^exp ( - -sin^/2) = —{[dtgr,f)) 
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We are reduced to the computation of {{dtgro}^)), and one gets 

mgrof)) = IcHl , 

with the factor 7c — 0.05 for our fractal model. This numerical estimate was 
already done by Carati et al. ( |2008j ). in providing the estimate (11) of page 
four. So also the fifth line of (llip is proven. 

We finally come to the last relation ffTTjl . One has (again taking only the 
largest terms) 

^ ^ qO) - X 

j 



H — - > Ma — ; COS zt^ 

(A 3 |qO) _ x| 

. «M ^(kfzS^ , M,ipi cos 20 . 
c4 ^ qO)-x ^qO)-x 

i 

which shows, in a very simple way, that one has {dtT)) = 0. For what concerns 
the variance of Pi the usual correlation argument shows that 

7rr 

As before, in the expression of {{dtT>Y) one can put, for example, 6' = 0, so 
that one gets 



jk 

] Ro , 



where Mtot = NM is the total mass of the galaxies, Rq = c/Hq is the radius 
of the universe, 70 a numerical factor which depends on the distribution of 
the galaxies. For our fractal model, a numerical estimate provides 7^ — 0.08. 
Now, using for Mfot the value given by (fT5|) . one eventually finds 



{(dtVY) = ^-fHl , 

which, inserted in the expression for cr|,^ , gives the second of the last line of 
ffTTj) . This completes the proof. 
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B The Wiener— Khinchin theorem 

We deal here with the problem of how one computes the autocorrelation 
(/(x)/(0)) if the autocorrelation (9x/(x)c}x/(0)) of the derivatives is known 
(see (19]) . We will show that if one assumes 



<ax/(x)ax/(0)> = ((ax/(0))Oexp(-Y) , (18) 

then one has 

(/W/(0)> = 1^ (l - ^""^'^'(1 + ^)) ((^x/(0))^> . (19) 

This gives formule (fTOj) . if we take for /(x) the function (yf^^. 

To show this, we make use of the Wiener-Khinchin formula (see Wiener 
[1930]; Khinchin [T9M] ) 

(/(x)/(0))= / |/(k)|V^-dk, (20) 

where /(k) is the Fourier transform of /(x). Now, the Fourier transform of 
9x/(x) is simply given by k/(k), so that, once (9x/(x)9x/(0)) is known, the 
spectrum |/(k)p (and thus the correlation) is easily found. 

Indeed, inverting the Wiener-Khinchin formula, the spectrum is given by 

k'\m\' = 77^, [ (9x/(x)9x/(0)>e-^'^-dx 

_ {{ ^■^^ ^) ) f exp f — exp(— i/cr cos 6'))r^ sin 6* dr d6' d0 , 

where, in the second line, we used polar spherical coordinates (r, 0, 6) in M.^ 
with the z-axis parallel to k. The integral is elementary, and one obtains 



27r2 (1 + /2F)2 
which gives the spectrum 

f 

2TT''k'' (1 + Pk''Y 



Now the correlation follows from the Wiener-Khinchin formula quite easily 
because, using spherical coordinates in the momentum space k with the 
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z-axis parallel to x, and the expression just found for the spectrum, 
relation fl2U]) gives 



n U k (1 + /2A;2)2 

The integral in the second line can be computed for example by using the 
method of the residues, thus obtaining relation f lT^ . 
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